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We show that it is possible to find maximal violations of the CHSH-Bell inequality using only 
position measurements on a pair of entangled non-relativistic free particles. The device settings 
required in the CHSH inequality are done by choosing one of two times at which position is measured. 
For different assignments of the "+" outcome to positions, namely to an interval, to a half line, or 
to a periodic set, we determine violations of the inequalities, and states where they are attained. 
These results have consequences for the hidden variable theories of Bohm and Nelson, in which 
the two-time correlations between distant particle trajectories have a joint distribution, and hence 
cannot violate any Bell inequality. 



I. INTRODUCTION 



It is well-known that the position operators of a particle at different times do not in general commute. This is the 
reason why the notion of trajectories cannot be applied to quantum particles. But non-commutativity is also a useful 
feature in some experiments. In particular, it is essential in experiments aiming at violations of Bell inequalities. 
In this paper we show that the non-commutativity of positions at different times is sufficient for getting a maximal 
violation of the CHSH-Bell inequality, and find the states required for this. 

Our first motivation for investigating this was the possibility of using such Bell experiments as a refutation of 
Hidden Variable Theories which do assign a distribution of trajectories to every quantum state: In such a theory the 
positions at all times have a joint distribution, and therefore cannot violate a Bell inequality. Hence their predictions 
must be in conflict with quantum mechanics and, most likely, with experiment. After completion of our work we 
found that this line of reasoning had already been followed by Correggi and Morchio |7] . We nevertheless include our 
discussion, and emphasize some additional points. Technically, the Bell violations found in [7] are for particles in an 
external potential, whereas we look at two free particles. 

Our second motivation for the present work is the endeavour of finding a loophole free Bell test based on homodyne 
detection in quantum optics. In such a Bell measurement, each detection must be a function of just a single field 
quadrature, which is mathematically the same problem as using functions of a single position variable. This is 
impossible with Gaussian states, because the Wigner function then provides a joint distribution. But with the new 
abundance of non-Gaussian states recently realized in the lab [H [22] there is a chance to find a feasible setup. Here 
the knowledge of the maximally violating states may be helpful, although only as a rough indication where to look. 
It would be even better to be able to start from a given state, and to identify the quadrature measurements giving 
the best violation. 

Our paper and our results are organized as follows: In Sect. |H]we briefly describe how our result contributes to 
the debate about hidden variable theories. In Sect. |III| we provide some general background concerning violations of 
the CHSH inequality. Here we introduce techniques related to the universal C*-algebra generated by two projections. 
These techniques are well known in the operator algebra literature, but as far as we know they have not been applied 
to simplify the theory of Bell inequality violations. In Sect. [IV] we outline how to get maximal CHSH violations from 
position measurements. There are three different settings: (1) We choose the "+"-event of each measurement as a 
position outcome in some finite interval. When for both Alice and Bob, di,d,2 are the interval lengths used for the 
first and second setting, m is the particle mass, and t the time separation, then the attainable violation depends only 
on the dimensionless parameter u = md\d,2 / '(AtK) . We show that maximal violation is attained for infinitely many 
values of u and also in the limit u — *■ oo. (2) When the "+"-event means that the particle is on the positive half-line, 
maximal violation can be almost achieved, up to an arbitrarily small error. Hence there arc singular states (i.e., states 
not given by a density operator) for which maximal violation is attained. These are necessarily dilation invariant, up 
to a quadratic phase. Finally, (3) we look at periodic sets. It is well-known [5], that periodic functions of position and 
momentum commute, if the product of the periods is 27rfi,/integer. Translated to the setting of a free particle with 
time difference t between position measurements with period pi and p2, we find the two measurements to commute 
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whenever u^ 1 £ Z, where u — mp\p2 / (2itth) . Of course, in that case no violation of a Bell inequality is possible. 
However, we show that this situation is very unstable, i.e., that the maximal violation jumps from zero to a finite 
value for u arbitrarily close to an integer. In an appendix we have collected some technicalities regarding case (2). 

II. THE BOHM-NELSON THEORY 

When it first appeared, Bohm's hidden variable extension of quantum mechanics [5] met much opposition from the 
mainstream physical community because it appeared to violate some basic tenets of quantum theory. Heisenberg, in 
his paper introducing the uncertainty relations, had convinced the physics community that the notion of trajectories of 
individual particles had no place in the theory. There was even a theorem by von Neumann showing hidden variables 
to be impossible. To its proponents, Bohm was seen to restore Realism to physics, giving a complete moment-to- 
moment account of where all the particles of a complicated quantum system really were. In part, these were also 
the motivations of Edward Nelson for creating his "stochastic mechanics" [14] . In addition, he endeavored to give a 
derivation of the theory, which was at the same time a derivation of quantum mechanics itself. Both theories are 
embedded in a family of such theories parameterized by the diffusion constant in units of h [9] , with Bohm's theory 
corresponding to and Nelson's to 1. In the limiting case of infinite diffusion constant, we find a theory in which 
positions at different times are just taken to be independent |25j . For the conceptual problems we discuss here, the only 
salient feature of all these theories is that they provide a joint distribution for all particle positions at all times, such 
that the equal time probabilities for particle configurations exactly reproduce the quantum mechanical probability 
distribution IV^ 2 - This is the basis for the claim that Bohm's theory is empirically equivalent to quantum theory. At 
the very least, this agreement reassures us that some aspect of these "real" trajectories is correct. 

On exactly the same footing, let us look at another quantity, which makes sense quantum and hidden variable 
theories alike, namely the two-time correlations between distant non-interacting, but possibly entangled, subsystems. 
Of course, in the quantum case the positions of the same particle at different times do not commute, so quantum 
mechanics has no joint probability for these. But for correlations between different particles there is no such constraint, 
and we can directly compare the quantum predictions with the two-time correlation functions from the Bohm-Nelson 
theory. As we show below, the quantum and the Bohm-Nelson predictions turn out to be quite different (they also 
disagree between the Bohmian and the Nelsonian variants). So if we take the agreement of one-time correlations 
with quantum theory as evidence that there is something right about these trajectories, we are now forced to admit 
that there is also something wrong with them. Certainly, this disagreement completely invalidates the argument of 
"empirical equivalence" between Nelson-Bohm theory and quantum mechanics. We could even stage an experimentum 
crucis on the basis of the explicit states and observables computed below. There is little doubt how these would turn 
out, probably not even for the staunch defenders of these theories. So our argument in some sense refutes the 
Bohm-Nelson theory. 

Of course, we are aware that the Bohmians and Nelsonians know about this disagreement, and will not be impressed 
|26j . The simplest position is to include the collapse of the wave function into the theory (2j [16]. Then the first 
measurement instantaneously collapses the wave function. So if agreement with quantum mechanics is to be kept, the 
probability distribution changes suddenly. There is no way to fit this with continuous trajectories: When the guiding 
field collapses, the particles must jump. While the glaring non-locality of this process may be seen as just another 
instance of implicate order, it introduces an element of unexplained randomness, and demotes the Bohm equation (or 
Nelson's Fokker-Planck equation) from its role as the fundamental dynamical equation for position. 

This may be the reasons why many Bohmians adopt a strongly contextual view. In this view one has to describe the 
measurement devices explicitly in the same theory, so all trajectories depend on the entire experimental arrangement. 
Therefore the trajectory probabilities in two experiments, in which the measurements on particle A happen at different 
times, have no relation to each other, not even for trajectories of particle B. So the two-time correlations computed 
from the 2-particle ensemble of trajectories are never observed anyhow, and hence pose no threat to the theory. The 
downside of this argument is that it also applies to single time measurements, i.e., the agreement between Bohm- 
Nelson configurational probabilities and quantum ones is equally irrelevant. The naive version of Bohmian theory 
holds "position" to be special, even "real" , while all other measurement outcomes can only be described indirectly 
by including the measurement devices. Saving the Nelson-Bohm theory's failure regarding two-time two-particle 
correlations by going contextual also for position just means that the particle positions are declared unobservable 
according to the theory itself, hence truly hidden. 

In this consistently contextual version of the theory, there may still be those "real" trajectories, but they are only for 
the eyes of Bohm's Demon, or some such hypothetical creature. No physical interaction, not even an "ideal position 
measurement", will reveal them to the mere human. This certainly explains the apparent paradox that Bohmians on 
the one hand place so much value on being able to say where the particles really are, but are, on the other hand, so 
remarkably uninterested in actually computing trajectories. But when the interest in the real trajectories is gone, the 
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only gain from the whole theory seems to be a pro forma justification for saying that the hand of a voltmeter is really 
somewhere. The mountain in labor gave birth to a mouse. 

III. GENERAL STRUCTURE OF CHSH VIOLATIONS 

In this chapter we look at the general problem of finding the maximal quantum violations of the Bell-CHSH 
inequality, when the measurements of Alice and Bob are given. All this is well-known, but since we need it several 
times, it may be useful to state the criteria in a compact form. 

Each of the measurements in the CHSH setting is a POVM with outcomes +1,-1, which means that it is charac- 
terized by two positive operators F± with P + + P_ = I. We can parameterize such observables by the single operator 
A = F + — F_ , which gives the expectation of the outcome, and satisfies — 1 < A < I. Then F± = (1 ± A)/2. We will 
assume the measurement to be projection valued (i.e. F± — F±), which is equivalent to A = A* and A 2 = I. In the 
CHSH setting Alice chooses two such measurements, Ai,A 2 , and Bob chooses B\,B 2 . Since their respective labs are 
widely separated, they can make their choices independently, and we may take A and B as acting on the respective 
tensor factors of the Hilbert space Ha ® associated with the combined system. There is a quantum state p of this 
system, in which correlations Tr pAiBj can be determined. The CHSH- correlation is the linear combination of four 
such terms, which is the expectation of the operator 

T = A 1 (g>(B 1 + B 2 )+A 2 (g> (Bt-B 2 ). (1) 

Given the operators Ai,Bj, the supremum of the CHSH-correlations Tr pT attainable with quantum states p are 
given by the largest elements in the spectrum of T, and since we can invert the outcomes of Alice's measurements 
{A4 i— » — Ai) we are equally interested in the most negative expectations. To summarize, we are looking for the 
operator norm ||T||. For arbitrary operators X we have ||A|| 2 = ||A*A||, and since T — T* we have ||Tj| = ^/||T 2 ||. 
A simple algebraic computation using the properties of the operators A4, Bj stated above gives 

T 2 = 4(1 + A 3 ® S 3 ), (2) 

where we have denoted e.g. A 3 := (2i)~ 1 [Ai, A 2 }. Since A 2 = I, A\ is unitary, and we have -A1.A3.A1 = — A 3 . Hence, 
the spectrum of A3 is symmetrical around zero, with maximum equal to || -4.3 1| . This gives a compact expression for 
the maximal attainable CHSH-corrclation, namely 

||T|| = ^4(l + \\A 3 \\\\B 3 \\). (3) 

In particular, when either Alice's or Bob's measurements commute, so the norm of the corresponding commutator 
vanishes, we get ||T|| = 2, i.e., the CHSH inequality is satisfied. On the other hand, since || [A 1; A 2 ]\\ < 2 1 1 ] | |[A 2 | = 2, 
we have ||T|| < y8 = 2v2, which is known as Tsirelson's inequality. For our purposes, the main gain from (|3| is that 
the determination of the maximal violation is reduced to the estimates of commutators, which can be done separately 
for Alice and Bob. 



A. The algebra generated by two projections 

Note that both on Alice's side (and similarly on Bob's) only two projections Pi = (I + Ai)/2 are relevant. Let 
A(P\ 1 P 2 ) denote the algebra generated by the two projections Pi and P 2 together with the identity I. It turns out 
[U Uni HH HH] that this can be understood completely in terms of 2 x 2-matrices. Indeed, we observe that 

C= I -Pi -P 2 + Pi P 2 + P 2 Pi (4) 

satisfies 

CPi - PiP-Pi - PiC, (5) 

and a similar relation for P 2 . Therefore C commutes with the generating projections of the algebra, and hence with 
all of A(Pi,P 2 ). The central element C satisfies < C < I, because C = 1 - (P x - P 2 ) 2 , and 



-1 < -P a < (Pi - P 2 ) < Pi < 1, 
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so < (Pi — P2) 2 < I. Clearly, C — means that Pi and P2 are orthogonal, whereas C = 1 means that Pi and P2 
are equal. At these extremes, [Pi,P2] = 0. More generally, the commutator satisfies 

[P 1 ,P 2 ]*[P 1 ,P a ] = ~[P 1 ,P 2 } 2 = C(1-C). (6) 

Hence the largest norm for the commutator square is 4, attained at C = |l, where indeed the operators A; = 2Pi — 1 
appearing in the CHSH-incquality also attain their maximal commutator norm 

||A 3 || = HP^H = 2||[P 1 ,P 2 ]|| = 2V||C(1-C)||=1. (7) 

Now we can express every element of A(P\, P2) as a linear combination of the four terms 1, P\,P 2 , P\P 2 , multiplied 
by suitable polynomials in C. It is convenient to choose another basis, in which multiplication becomes ordinary 
matrix multiplication, and such an isomorphism can be implemented at the Hilbert space level (see e.g. (TUl [H]): 
First we split off the null space of all commutators, i.e., H.q = {</>|C(J — C)<p = 0} = keryl.3. On this space, which is 
clearly an invariant subspace of the two projections, all four eigenvalue combinations of two commuting projections 
are possible. Going to the orthogonal complement TCq, we put JC := {<fi G Ho\Pi<fi — </>} — PiHq, and let H denote 
the restriction of the central element C to this invariant subspace, i.e. H = C|jc = PiP^PlIk:- Then we define 



by 



V : Hi -> K <8> C 2 (8) 



Vcp = iPicj>® + —====PiP 2 {l - PO0® |-) x , (9) 



where li)-,^ = (±i)- C> ne can readily check that this map is unitary. Operators on JC ® C 2 can be considered as 
S(/C)-valued 2 x 2-matrices; this gives 

Ai| W J. <=±I®(Ti, 

A 2 \ n< L = a(H) (g> ai + /3(H) (g> a 2 , (10) 

where a, j3 : [0, 1] — > R are given by a(/i) = 2/i — 1, /3(/i) = 2-»//i(l — h). The point of this decomposition is, of course, 
that in the matrix entries we only have functions of the central element C or rather its compression H . 



B. Attained maximal violations 



We now use the detailed form ( 10 1 of the operators to get better information about the states where large violation 
is attained. Both algebras now have a central element, C a and Cb, respectively, giving the compressions Ha and Hb- 
When these are fixed ft,- numbers, the four operators Ai,Bj are completely fixed 2 x 2-matrices, and we can explicitly 
find a state on C 2 <g> C 2 maximizing the violation. In general we can do this maximization at every pair of values, 
which by ^ gives the function 



p(h A ,h B ) - 2^1 + VMi -^a) VMi - ^b) (11) 

plotted in Fig. [T] Given the joint probability distribution of Ha and Hb, the largest attainable CHSH-correlation 
will be the expectation of (111 with respect to this distribution. Obviously, for large correlation we want to choose a 
joint distribution which is concentrated as near the point Ha = hs = \ as possible. 
This leads to the following three cases: 

1. When i is an eigenvalue of both Ha and Hb, we will chose a maximizing vector from these eigenspaces. Then 
the CHSH-correlation will be exactly 2\[2. 

2. When 1 lies in both spectra, but for one of these operators is not an eigenvalue (i.e., lies in the continuous 
spectrum), the CHSH-correlation for any state represented by a density operator in the given Hilbert space will 
be strictly less than 2^/2, but can be chosen arbitrarily close to this value. 

3. Finally if \ is not in the spectrum of H,\ or Hb- the CHSH-correlation is less than 2\f2 — e for some s > and 
all states. 



FIG. 1: Maximal CHSH-correlation as a function of the central parameters of Alice and Bob 



To characterize the structure of the maximally violating states, we would now like to extract from the case 1 as much 
information as possible about further expectation values, including those not directly measured in a Bell experiment. 
Similarly, in case 2 we are interested in the limits of expectation values Tr(p n A), for p n a sequence of density operators 
with asymptotically maximal violation. It is convenient to treat these two cases on the same footing by choosing a 
convergent subsequence of the p n in the weak*-topology, and thereby find an exactly maximally violating limiting 
state. This is possible if we extend the notion of "states" from density operators to arbitrary expectation value 
functionals lo : B{TLa ® 7~Lb) —> C Of these we only require linearity, positivity and normalization, so they are states 
in the sense of C*-algebra theory. The state space of a C*-algebra is compact with respect to "convergence of all 
expectation values", so convergent subsequences in this wider setting always exist. Of course, in case 2 a sequence p n 
can converge only to a singular state and not a "normal" one, given by a density operator. The singular state is never 
unique, because fixing such a state is the non-commutative analog of explicitly defining a free ultrafiltcr. However, as 
will be seen below, all these states may well agree on some observables of interest. 

Geometrically, the CHSH expression with fixed Aj, Bj is an affine functional on the state space of B(TCa ® Wb). It 
reaches its maximum at an extreme point, so it is not surprising that this entails some special relations. The typical 
tool here is the Cauchy-Schwarz inequality in the form |w(X*F)| 2 < u}(X*X)cu(Y*Y). Hence if we know that the 
expectation of a positive operator like X*X vanishes we can conclude that oj(X*Y) = lo(Y*X) = for all Y. This 
approach has been applied [21 to the CHSH inequality by writing 2 \/2 — T as the sum of several operators of the 



form X*X. Here we can achieve similar results by looking at the explicit form (10) of the operators A i: Bj after the 
transformation ([£]). 

Under this transformation Alice's Hilbert space becomes H A © (ICa <£> C 2 ), so that Ha is an operator on K-a, and 



the given operators A\,Ai take the form (10 1. Of course, an analogous decomposition holds for Bob. The projections 
onto these subspaces as well as the spectral families of the operators Ha,Hb commute with all A i: Bj. Suppose we 
take a joint spectral projection P £ of the commuting operators Ha,Hb corresponding to a set of distance e to the 
point (i, |). Then TP e is strictly smaller than 2-^/2, so a maximizing state must have uj(P e ) = 0. Hence a maximally 
violating state must vanish on Hq^Hb and Ha^TLq, and its restriction to B(ICa® ICb) must be a state Si i giving 
the sharp values | to both Ha and Hb, in the sense that 

5 U ((H A -lir)=S u ((H B - 1 -ir)=0. (12) 



At such a point we can set a(h) = and (3(h) = 1 in formula ( 10 1, and its analogue for B\, B^ and just consider the 
maximization of the CHSH expression with fixed qubit operators A, = B t = a^. Since the maximum of T for these 
operators is attained at the unique pure state 

*o = 4=( e_W4 l ++> + |- -))eC 2 ®C 2 = C 4 (13) 
v2 

(where |±) are the eigenvectors of 03), we conclude that the overall state must be of the form 

u = <g> |*o)(#o|. (14) 

It is clear, that conversely, any state of this description will be maximally violating. 

In case 1, the explicit form for the maximally violating pure states in the original representation is now 

* = -±-[e-™l 4 e A <+ <g> e B '+ + e A >~ ® e s <-], (15) 
v2 
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where e.g. 

e A >± := VZUa ® |±)) = {-iPf ± (1 - P 1 A ))(V2P A f), (16) 



and /a is a normalized eigenvector of Ha belonging to the eigenvalue |. In case 2, we choose a sequence (/^) of unit 
vectors such that \\HAf A — hf A \\ ~ > (usually called approximate eigenvectors). Defining ty n using f A and f£ as 



in (15 1, we get the asymptotic maximal violation hrn n _ >00 (\I r 7l |T'\I r ri ) = 2y2; the states |\E , n )(\E' n | approximate some 



singular state of the form (14 1. This systematic construction of (approximate) maximally violating wave functions 
will be used in the next section. 

An interesting corollary of the above structure is the cryptographic security of maximal CHSH correlations. We are 
then interested in the possible correlations between the observed data and the measurements made by an eavesdropper 
"Eve" in a separate lab. The measurement of Eve is then described by an operator E commuting with all the operators 



Ai, A 2 , Bi,B 2 used by Alice and Bob. Hence E lives on the tensor factor Ka®^Bi and from the form ( 14 1 of the state, 



it is clear that Eve's results are independent of Alice's and Bob's. This is summarized in the following proposition. 

Proposition 1. Let lu be a state maximally violating the CHSH inequality on operators A\, A 2 , B\, B 2 . Let p be 
a non- commutative polynomial in four variables, and set P — p{A\, A 2 , Bi, B 2 ) and P° = p{A\, A 2 , B®, B 2 ) , where 
A® = ai ® I and Bf = I ® <r,. Then for any operator E commuting with all Ai, By. 

co(EP)=u J (E)(^ \P ^o)- (17) 



IV. POSITION MEASUREMENTS AT DIFFERENT TIMES 



We now proceed to the case of position measurements. The Heisenberg picture position operator of a massive, 
freely evolving nonrelativistic particle with mass m is given by 

Qt = Pt/m + Q, t e R, (18) 

where Q and P are the standard position and momentum operators, acting in the Hilbert space L 2 (JBL) (in particular, 
P = —iH-j^). Concerning measurements of Qt (position at time t), we are only interested in recording whether the 
outcome lies in a fixed interval A C 1, in which case we assign the value "+1" to it; otherwise we label it "—1". The 
corresponding two- valued observable is 2xa(Q*) — 1- 

We consider the case where Alice makes measurements with one particle, and Bob with another one; let A\ := 
2xAi (Q A ) — 1 and A 2 := 2xa 2 (Q A ) — 1 be the position measurements for Alice's particle at time zero and time t > 0, 
with intervals Ai and A2, respectively, and let Bi be the similar ones for Bob. For simplicity, we suppose that both 
use the same measurement intervals, same time t, and particles of same mass m. Now we are in a situation discussed 
in the preceding section. 

Since the operators are identical for both parties, we consider only Alice's part and drop the associated index when 
there is no confusion. The two projections are now Pi — XAi(Q) and P 2 = XA 2 (Qt)- We begin with the fact that the 
pair (Q,t -1 mQ t ) is canonically conjugated, and therefore unitarily equivalent to the pair (Q, P), the unitary operator 
in question being simply e^ 2 *^ lm Q"/ fi . With this equivalence, Pi ~ XAi(Q), P2 — Xt- 1 mA 2 (P)> m the following, we 
will simply replace the with these operators. 

The idea is to consider the possibility of maximal violation of the CHSH inequality for three types of concrete 
choices for the localization intervals Aj, exhibiting different commutativity behavior of the position and momentum 
projections Pi and P 2 [5]: 

(1) For compact intervals, Pi and P 2 are partially commutative, i.e. ker[Pi,P2] = ker^3 is neither {0} nor TL. 
Indeed, there are common 0-eigenvectors of Pi and P 2 . 

(2) For half-lines, the projections are totally noncommutative, i.e. ker[Pi,P2] = {0}, and K. = P 2 (Ai). 

(3) For periodic sets, the periods can be chosen so that Pi and P2 are commutative, i.e. ker[Pi,P2] = TL. Then 
K = {0}. 

(The full characterization of commuting functions of Q and P is given in 6 ; for a generalization to Abelian groups, see 
[24].') The projections apparently depend on various parameters Ai, A2, t, m; however, as the dilations are represented 
by unitary operators, the only relevant parameter is the scale of the Q-interval relative to the P-interval. In case (2) 
there is no specific scale, because the projections are invariant under dilations; hence the structure of the Bell inequality 
violations does not depend at all on the parameters. In cases (1) and (3), Q- and P-sets are characterized by lengths 
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l\ and mf _1 ^2, respectively, where the k are proportional to the lengths (case (1)) or periods (case (2)) of the sets Aj. 
If we fix the unit of position as l\ (thereby making the position variable dimensionless) , the unit of momentum will 
be fill 1 ] m t nese units, the above characteristic lengths are 1 and 

u = ^h~> (19) 

respectively. We can equally well fix the unit of momentum as mi -1 ^, in which case the unit of position is tft/(mZa); 
the characteristic lengths are then u and 1, respectively. Hence the only relevant parameter is the dimensionless scale 
u. For technical reasons, we will use the first choice of units in case (1) and the second in (2). For both choices of 
units, the associated operators are dimensionless; we will denote these by Q and P. 



A. Compact intervals: partially commutative case 

Here we let Aj ^ be a compact interval for i — 1,2. As already mentioned, ker^ is nontrivial; however, 
P\(7~t) = L 2 (Ai) c Hq (see e.g. [5]), so the relevant subspace /C is just L 2 (Ai). It is convenient to choose the length 
scales as k := di/2, with di the length of A*; passing to the units where l\ is 1 as discussed above, we see that the 
relevant operator H is unitarily equivalent to 

H u := X[-i,i](Q)X[-u,u]( P )X{-i,i](Q) 6B(i 2 ([-U])), 



where u is given by (19 1, i.e. u — md\d2 / '(4th) . (This equivalence can be seen easily by first applying the usual 
translation and "velocity boost" unitaries with appropriate shifts to center the intervals to the origin, and then 
dilating by di/2.) 

The structure of H u has been extensively studied because of its relevance in band- and timelimiting of signals (see, 
for instance [8j pp. 21-23], [TTJ pp. 121-132], [20], or the original papers by Landau, Pollack and Slepian [T3l [T7l [20] 1. 
We summarize the relevant mathematical facts briefly in the following paragraph. 

The operator H u is explicitly given by 

(H u <p)(v) = f M<v-w)) y e L*([-l,l]), 
J-i n(v-w) 

from which it follows that H u commutes with the differential operator [(1 — w 2 )^;] — u 2 v 2 that determines the 
angular part of the wave equation in prolate spheroidal coordinates. This differential operator has a complete or- 
thonormal set of eigenfunctions ipn ^ L 2 ([— 1, 1]), n = 0, 1, . . ., called angular prolate spheroidal wave functions. In 

the notation of [19] we have ipn( v ) = \J n + ^Psn(f, u). The corresponding eigenvalues X n (u) of H u are 

X n (u) =2utt- 1 S^(1,u) 2 G (0,1), n = 0,1,2,..., (20) 

where Sn\',u) is the radial prolate spheroidal wave function of the first kind. In particular, X n (u) depends continu- 
ously on u. In addition, we have 1 > X n (u) > X n +i(u) > for all n and u. 

Now | G cr(H) exactly when u is chosen so as to make X n (u) = | for some n. Since lim^-xx, X n (u) = 1, and 
lim^^oo X n (u) = for fixed n (see |20j). it follows by continuity that for each n we get at least one value u n G (0, 1) 
with X n (u n ) — \. On the other hand, H u < H u i if u < u' , so each X n (u) is an increasing function of u, and u n is thus 
uniquely determined. Since for given n, we have X n (u) > A, i+1 (u) for all u it follows by continuity that u n < « n+ll 
i.e. the sequence (u n ) is increasing. 

Figures [2] and |3] show the u-dependence of the largest eigenvalues, as well as the relevant commutator norm H^H, 
obtained from the above representation. The critical values u n can be computed numerically; the smallest two are 
approximately u Q sa 0.849 and U\ ~ 2.381. 

Summarizing, a state that maximizes the violation of CHSH inequality for position measurements exists if and only 



if both Alice and Bob adjust their parameters in such a way that ( 19 1 holds with u one of the critical values, say u n 
for both Alice and Bob [27]. Using ( 15 ), the corresponding wave function ^ G L 2 (R 2 , dqA, dqs) can then be expressed 



quite explicitly. For this we need the functions e ± of (16); but now v^^V-'n" = \J \( n + \)^ s n(', u n) G L 2 (R), so 
e ± (q) is simply this spherical function 
function is discontinuous at the lin 
origin, i.e. Aj = [-di/2, di/2], we get 



e ± (g) is simply this spherical function, multiplied with — i if q G Ai, and with ±1 otherwise. In particular, the wave 
function ^ is discontinuous at the lines qA — ±1 and qg = ±1. In the case where the intervals are centered at the 
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FIG. 3: The norm of the operator A3 as a function of u 



y(q A ,q B ) = C (q A ,q B ) -J-(n + l)Ps n (2q A /d u U n )Ps n (2q B / d u u n )e ie ^ B \ (21) 



where 



. v 'l + l/\/2, q A ,q B S A 1)0 r q A ,q B $ A t 

Lo{qA,qB) = \ 1 — (22) 

1 — 1 / v2, otherwise 



is real, and the phase Q(q A , q B ) is given by 



Q(qA,q B ) = -2u n (q 2 A + q%)/(d 1 d 2 ) + { q A ,q B £^ (23) 

bQ+ir/2, otherwise. 



0± = =Farctan((%/2± 1) _1 ). 
Figure [4] shows the picture of the simplest choice for the wave function. 



B. Half-lines: totally noncommutative case 



Here we set A, = [xi, 00), i = 1,2, i, 6 K; then ker ^3 = {0}, and we have the tensor product representation for the 
full operators. As before, JC — L 2 (Ai). We first apply the unitary shifts that transform the situation to the dilation 
invariant case A, = [Q, 00), i = 1,2; then the units of position and momentum play no role, and A\ = Sign(Q) and 
A2 = Sign(P). In particular, the spectrum of H does not depend on any of the parameters. 
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FIG. 4: The maximally violating wavefunction ( |21[ ) with n — (u = uq w 0.849), Ai = A2 = [— 1, 1]. The plotted function is 
^(qaj^b) without the complex phase factor e ie ^ ,j4,9B '; the phase 0(gA, </b) is shown in color, going through red, yellow, green, 
cyan, blue, magenta, and again red, as the value varies from to 2n. Note the discontinuity lines qa = ±1, <1b = ±1 of the 
factor Co(qA, Qb) (see equation (22 1), marking the measurement interval. 



The transformation V : L 2 (R) — > L 2 ([0, 00)) ® C 2 given by (|9) has now a particularly simple form, namely 

V r ^ = t[x [ o,Do)(Q)0® |+)i + nx(_oo,o](Q)0® Hi]. ( 24 ) 

where II : L 2 ((— 00, 0]) — » L 2 ([0, 00)) is the parity operator. However, we still have to determine the spectrum of the 
operator H — X[o,oo)(Q)X[o,oo)(-f )X[o,oo)(Q)> acting on K — L 2 ([0,oo)). Here it is convenient to utilize the dilation 
invariance of the projections; we seek a unitary operator W that diagonalizes the dilation generator D := QP + PQ 
by way of 

WDW* =Q®I c i (25) 

Such a unitary is obtained by first separating the positive and negative half-axis by using the V above, then expanding 
X 2 ([0,oo)) into the full L 2 (R) via the unitary operator U+ : L 2 ([0,oo)) -> L 2 (R), where 

(U +V )(\) = V2e x V (e 2X ) 1 (26) 

and then applying the Fourier-Plancherel operator F. Indeed, the operator (L/+<B> Ic 2 )V transforms D into P(8>lc 2 , so 
we get ( 25 1 by setting W := —i(FU+ ® Ic 2 )V '■ L 2 (R) —> L 2 (M.) (g> C 2 (where the factor —i is chosen for convenience.) 
This unitary operator approximately diagonalizes Pi and P2 simultaneously, meaning that we get explicit form for 
the functions of H appearing in ( JToj ) (see Appendix). We can then explicitly compute H = P1P2P1: this gives 

H= -(J + tanh(|7rQ)), (27) 

acting on the space L 2 (R) = U + L 2 ([0, 00)). 

It follows from (27 1 that the spectrum of H is purely absolutely continuous, and contains the point g. Hence, 
maximally violating states lu exist, are of the form (14 1, and each of them is singular. Note that the tensor product 
representation space is now L 2 (M. 2 ) ® C 4 . 
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It is possible to further specify the properties of the restriction of u> to the first tensor factor L 2 (M. 2 ). According 
to (27 1, spectral projections of H associated with intervals around \ correspond bijectively to those of Q around 



b2 



0. Hence, equations (12) imply that 8i f i {f{Q A , Qb)) — /(0,0) for any bounded measurable function / 
continuous at (0,0). With this information, we can now go back to the representation where T is given by ([lj with 
A\ = Sign(QA), A2 = Sign(P/i), and Bi similarly; using (251, as well as the Schwarz inequality, we see that any 
maximally violating state u> satisfies 

w(/(D A , D B ) X) = u{X f(D A , D fl )) - /(0, 0)u(X), (28) 

if / : M. 2 — > C is a bounded measurable function continuous at (0,0), and X G B(L 2 (M. 2 )) is arbitrary. In particular, 
each maximally violating state is invariant under dilations in this representation. Since Q and P transform covariant 
under dilations, this means that maximally violating states are concentrated on (0, 0) and infinity, in both position 
and momentum representations; the precise statement is the following observation. 

Proposition 2. Let u> £ B(L 2 (JSL 2 ))* be a dilation invariant state, and let f : R 2 — > C be a continuous function 
vanishing at the origin and infinity. Then 

w(/(Qx,Qb))=w(/(Pa,Pb)) = 0. 

Proof. Set K(ri, r 2 ) := B(r 2 )\B(ri), where B(r) is the open ball in M 2 of radius r centered at the origin. By dilation 
invariance, w(xs(r)(QA, Qb)) = v{Xb(i)(Qa,Qb)) for all r > 0, so for r 1 < r 2 we get u)(xK( ri ,r 2 ){QA, Qb)) = 0. From 
the positivity of oj it follows that Lj(f(Q A , Qb)) — for any bounded measurable function / : K 2 — > C supported in 
K{r\, r 2 ). Now if / : K 2 — * C is continuous and vanishes at both zero and infinity, then linin^oo \\xK(i/n.n)f~ f\\oc = 0, 
and hence u>(f(Q A , Qb)) = by the norm continuity of w. The case with is similar. □ 

We now wish to find wave functions approximating the maximally violating singular states lu. Since any such 
state is dilation invariant, one can expect that the approximating wave functions in L 2 (M. 2 ) would basically look like 
l/^y\xy\, but with some regularization at the coordinate axis and infinity. (Here \j\fx comes from formally solving 
the "eigenvalue equation" Dip — 0.) 

In order to construct such approximating wave functions, we proceed as described in the preceding section. The 
approximate eigenvectors of ( p7| ), corresponding to the point | £ a(H) are of the form g e — (2e)~ ? g(x / (2e)) , 
where g E L 2 (M) is an arbitrary unit vector, and e > is small. Hence, the corresponding vectors for the original 
H = X[o : oo)(Q)X[o,oo)(-P)X[o,oo)(<2), acting on L 2 ([0,oo)), are / e := iU+F* g e ; explicitly, they are of the form 

f e (x) = iJ^f(elnx), x > 0, 



where / G L 2 (M) is an arbitrary unit vector. Hence, the wave functions we are seeking are given by (15 1, with 
ef — V*(f e ® |±)) (for both Alice and Bob). These can now be obtained from ( 24 1 : 



<{x) = ^M\x\), e~{x) = -^Sign(x)/ e (M). 

Hence, 

* e (qA,q B ) = -^(e-^ 4 + Sign( qA qB)) ^=f(eln\q A \)f(eln\q B \). (29) 
2y2 Vl9A9s| 

This same formula appears in [lj; however, the paper does not seem to contain any systematic derivation for the 
result. 

The approximating wave function in the original representation, where both A4 are position measurements, is then 



0?a,9b) ^ ^ e {q A ,q B )e ^' 



^ 1 (q 2 A +q 2 B ) 



where the measurement intervals are [0,oo) for both time zero and t. 

We close this subsection by demonstrating that the singular states that can be approximated by the wave functions 



( 29 1 actually depend on the regularizing function /, even though they are all maximally violating and dilation invariant. 
One property of a dilation invariant state u that we can easily determine is the expectation value u>q := uj(h(QA, Qb)), 
where h : K 2 — > M is bounded, measurable, and continuous at the origin, with /i(0,0) = 1. By dilation invariance, this 
does not depend on h, and describes the "weight" of the state at the origin in the position spectrum. Note that by 
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Proposition [2] each dilation invariant state is concentrated at zero and infinity. However, the distribution of weight 
between these points is not fixed: by direct calculation using (29 1, we get 



W = l ^(^e\X[-a,a]x[~a,a](QA,QB)^e) = Hm 

c_j.ll L ' i ,\\ 



\f(x)\ 2 dx 



\f(x)\ 2 dx 



(30) 



which may attain any value in [0, 1], depending on where / is concentrated. 



C. Periodic sets: commutative case 



Here we simply want to make a remark about the commutativity, without trying to analyze the periodic set case 
systematically. Consider sets of the form 

A i=Pi Z + [Q,Pi/2], 

where Pi > are the periods. Now we choose the length scales as := pi and l 2 := p 2 /(2n). Passing to the units 
where mt^ 1 ^ is 1 (see the beginning of the section) via the associated dilation, we get A\ _ g u (Q), A 2 — g2^(P), 
with 



9v{q) 



1, gew([0,|)+Z), 
— 1 , otherwise 



and u again given by (19 1; u — mp\p 2 / '(2irtK) . The reason for the choice of units is that A\ and A 2 (or, equivalcntly, 
Pi and P 2 ) commute if u^ 1 is an integer. This can easily be seen by noting that g v (x) = Sign(sin(27ra/t;)), and using 
the commutation relation for the Weyl operators. Moreover, the converse is also true; see the general characterization 
of commuting functions of Q and P [5J [M] • In the commuting case the spectrum of H contains only the points 
and 1, and we have A3 = B3 = 0. The CHSH-inequality is then actually satisfied for all states, and the situation is 
classical. 

It is interesting to observe that when the parameter u^ 1 is slightly perturbed from an integer, the commutator 
norm \\A^\\ discontinuously jumps to a nonzero value, which is large enough to allow a violation of Bell's inequality 
In order to show this, we take u = 1 + e, with e > 0. Now g 2v has the Fourier expansion 

S2„(p) = $>n^Le^, 

where c_„ = — c„, and J2 neZ \c n (u)\ 2 = \s u (p)\ 2 dp = 2tt. For each < e < 1/2 choose a unit vector tjj e S L 2 (M) 
with support in [0, e). Then g\ +t (Q)ip t — ip e , and 

~ri v 27r 

where the series converges in L 2 (M.) because the terms are orthogonal. Since (e mP ip e )(x) — ip e (x + n), we get 

\-e m ^ip c , -|(1 + 1/e) < n < 

It follows that 

[g2n(P),Sl +e (Q)]^e = 2 ]T ^C„e mP ^ + ^, 

-l/(2e)<n<-l * 

where ip t is orthogonal to the sum. Hence, 

i[ 52 .(p), Sl+e (Q)]ii 2 > - J2 i c «i 2 - 

l<n<l/(2e) 

Here the right hand side tends to 2 as e — > 0, which means that \[g 27r (P) , gsi +e (Q)]\\ > y2 for all sufficiently small 
e > 0. Hence, | j ^4.3 1 1 jumps discontinuously from to some value larger than 1/V2, corresponding to the Bell correlation 
V6. 
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Appendix: Almost simultaneous diagonalization of Sign(Q) and Sign(P). 

The diagonalization is provided by the unitary operator W : L 2 (R) — > L 2 (R) ® C 2 defined in Section 
can be written in the form 



IV B 



This 



w<i> = (fu+ ® s)( X[0 ,oo)(Q)<f> ® 1+) + nx ( _oo,o]0 ® I-)), 

where |±) are the eigenstates of (T 3 , 5 is the Hadamard matrix S = 4= (} and J7+ is defined in (26). 
For Sign(Q), we get 

WSign(Q)W* =7® (7i 

by trivial calculation. The form of WSign(P)W* is not so obvious, and the following computation actually describes 
how to find a suitable W. Let cp G L 2 (R) be a Schwartz space function, with <p(Q) = 0. The set of such functions is dense 
in L 2 (M). By using dominated convergence twice, and then Fubini's theorem (noting that f_ \ip(x)\/\x\ dx < oo) we 
get 

(ip\Sign(P)(p) = liiri V / <p(x)K' t 5 (x,y)ip(y)dxdy, 

ei,e2=±l =i,«2 



with 



1 

2^ 



1 



Jeicc — i(x — y) 5e 2 y + i(x - y) 



where M. 2 ± £9 denotes the appropriate quadrant. Each kernel K% t2 is invariant under dilations, i.e. aKf e2 (x,y) = 
e2 (a~ 1 x,a~ 1 y) for all a > 0. Using this we can transform them into convolution kernels, and then diagonalize 
using the Fourier transform. Indeed, put £2 (A) := 2 AT* £2 (eie A , £2e~ A ), AgK; this gives 



*£, ea (A) 



<5e A 



2iG £l;£2 (A) J e - A + 2iG £l , £2 (A) 



where G±± = ±sinhA, and G± T = ± cosh A. For S > 0, each If* £2 is both integrable and square integrable, so we 
can put K% i ea {r)) := J* R e'^K^ C2 (A) dA. Then we compute 



tp(x)K> eue2 (x,yMy)dxdy = / dA (l/ £l ^)(A)« !£2 * V^)(X) 



= J dti{FV ei <p)(v)Kt uea ( V ){Vt a <p){ri) = (FK^|<, £2 7V e2 ^}. 

It remains to take the limit S — > 0+. To this end, first note that since e.g. |(5e A — 2iG± i=F (A)| 1 < (coshA) -1 for 
any 5 > 0, and since A i— ► 1/ cosh A is integrable, it follows that .K± T are bounded uniformly for S, and if± =F ( ? 7) 
exists with K± ^(rj) — lim £ ^o+ K± ^(v) pointwise by dominated convergence. Hence, the corresponding bounded 
multiplication operators on L 2 (M., drf) converge in the strong operator topology, giving 

l\m + (FV e MKi T V e2V ) = (FV ei v\Kl T FV e2 <p). 
Since K± ^(A) = =F(i7rcosh A) -1 are even functions, we have 



cos(i]X) 
cosh A 



dX = ± 



cosh(?77r/2) 



For diagonal elements, K± ± = =F(i7rsinhA) 1 , and the corresponding Fourier integral does not exist. However, the 
limit K± ±(rj) := lim,5_>04- K± ±( r l) exists pointwise, because K± ± is an odd function; in fact, 



K ± ±(v) = ~2i / sin(r)X)Ki ± (X)dX — > ±- 



sin(?7A) , /??7i"s 

/ dX = ±tanh 4r ), 
sinhA v 2 ; ' 



(31) 



13 



as 5 — > 0+, the singularity at the origin being canceled by the factor sin(?7A). We can now use e.g. the bound 
\K±,±(v)\ < M\n\ < M(r] 2 + 1), where M = f J °° A(sinh A)- 1 d\ < oo, and the fact that HQF^^ = \\PV ± tp\\ < 
\\ip\\ + 2\\Q(p\\ < oo, to conclude that 

lim (FV e MKl,±FV e2V >) = {FV t MKl.±FV t2 ^. 

The coefficient matrix is thus K (rj) — tanh(r77r/2)<73 — sech(?y7r /2)<72- Finally, taking into account the Hadamard 
matrix S in the definition of W, we get the result 

WSign(P)W* = tanh(Q7r/2) <g> a t + sech(Q7r/2) <g> cr 2 - 
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There are further generalizations. In fact, it is easy to build Markov processes following the quantum evolution of any 
observable 123 , so one could directly discuss everything in terms of spin variables rather than positions. However, this 
generalization is clearly against the taste of Bohmians, and we will ignore it as they have. 

A notable exception is the founder of stochastic mechanics, who abandoned the theory, when he realized some of its 
unphysical non-local features [T5] 

In order to simplify the form of the maximally violating wave function, we have assumed in the beginning that Alice and 
Bob use the same parameter values. Clearly, they could just as well use different values, in which case Alice must have 
it a = u n for some n, and Bob can have ub — u n i for some other number n' . 



